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1. Introduction 

It is well known that chaotic dynamical systems may exhibit deterministic diffusion: 
among the simplest examples one may recall chains of 1-d maps P], area preserving 
maps on the cylinder |2] and infinite domain billiards, like the Lorentz gas with finite 
horizon [3]. In many cases the mean square displacement (the average being taken with 
respect to initial conditions) displays a linear asymptotic growth 



and one may try to devise techniques to compute the diffusion constant D. In 
particular periodic orbits expansions have been employed jl] to evaluate D for one 
dimensional maps [H] , finite horizon Lorentz gas [HI and area preserving sawtooth maps 
[7j. Deterministic transport may also exhibit an anomalous character, with a mean 
square displacement deviating from a linear behaviour. Typically this is associated to 
a weakening of chaotic properties, like marginal fixed points |H] , the presence of regular 
islands in the chaotic sea ^ , or the opening of ballistic channels for Lorentz gas models 
[10]. Such weakly chaotic regimes are generically quite hard to study, as they miss the 
fast memory loss that typically characterizes strongly chaotic systems. Besides transport 
properties this lack of full chaos deeply modifies other dynamical features: for instance 
escape processes may be slowed down and time correlation functions acquire power 
law tails (in view of Green-Kubo formula ^I] this is closely connected to anomalous 
diffusion). Global properties of the transport process are investigated through the 
asymptotic behaviour of generalized moments 



Gaussian processes yield i/(g) = q/2 for all q; recently deviations from a constant 
slope for z/(g) have been referred to as "strong" anomalous diffusion ^2] (the reader 
may consult this reference for a wide survey of systems falling in this cathegory). In 
this paper we will introduce a cycle expansions technique for computing the spectrum 
z/(g) {q/2 G N) and apply the formalism to a number of examples: different chains 
of one-dimensional intermittent maps, and infinite domain billiards (Lorentz gas with 
infinite horizon). 

The paper is organized as follows: in section 2, to make the paper self contained, 
we review how periodic orbit expansions may be employed to study transport processes, 
section 3 briefly explains how full chaos generically results in a gaussian i'{q), while, 
finally, section 4 is devoted to the analysis of weakly chaotic systems yielding a non 
trivial z/(g) (parts of the results of this section were anticipated in jTHj)- 

2. Deterministic transport 




(1) 




(2) 



Chaotic dynamical systems are characterized by sensitive dependence on initial 
conditions, so that a meaningful theory must concentrate on statistical properties, rather 
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than on analysis of single trajectories. A central role in such a context is played by the 
Perron-Frobenius operator C which acts as: 

[£0](x)= / dySifiy)-x)<Piy) (3) 
Jx 

(we are considering for simplicity dynamical evolution described by a discrete mapping 
/). This operator may be interpreted as the evolution on densities (dual picture with 
respect to Koopman dynamics on observables) 

Pn+iix) = C[pn]{x) = dy 5{f{y) - x)p„(?/) (4) 



Jx 

(in particular the invariant density is an eigenvector of C, with corresponding eigenvalue 
A = 1). Notice that the kernel is singular, so that investigation of spectral properties 
generally requires sophisticated techniques jT^: we also observe that this operator enjoys 
the semigroup property: 

[C"^^]{x) = [C--'[C'<P]] (x) (5) 

The study of transport properties is conveniently carried out in terms of the generating 
function: 

GtiP) = (e^("*-"»)) (6) 

when again the average is over a set of initial conditions. We also define a generalized 
transfer operator where the singular delta kernel in is weighted by an appropriate 
factor (that maintains the semigroup property): 

[C^<P]{x) = [ dy e^^f'^y^-y^dUiy) - x)<P{y) (7) 
Jx 

The key idea is that if we are able to derive an analytic expression for the generating 
function, we can then extract all transport moments: 



(8) 

/3=0 



The generating function may be expressed as 

cm = j da;[4pij(x) (9) 

where pm is the density of initial conditions: the asymptotic growth of Gt{j3) in this 
way is dominated by the leading eigenvalue of the generalized transfer operator, as in 
equilibrium statistical mechanics of lattice systems, namely 

Gt{f3) ~ X\f3) (10) 

and then 



813' 



1 



(11) 

/3=0 

The main problem thus consists in computing the leading eigenvalue (and its 
derivatives): before turning to this point we have to describe the particular class of 
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dynamical systems we want to study: generally we are interested in spatially extended 
periodic systems where the dynamics over the unbounded phase space is generated by 
replicas of a map on a finite fundamental domain (like the Lorentz gas on a square 
lattice, where the dynamics on a fundamental cell corresponds to the Sinai billiard). 
From now on we illustrate the case of a one- dimensional map on the real line, even 
though the method is by no means limited to this context. The key property of the 
dynamical system under investigation is represented by the symmetry properties 

f{-x) = - fix) (12) 
f{x + n) = n + f{x) (13) 

for any n E 1^. The property (fT^ guarantees the absence of a net drift while with (fT^ 
we impose that the map on the real line is obtained by lifting a circle map f{6) defined 
on the unit torus as 

/>) = /WLodi eeT=[0,l) (14) 

We may thus split the / evolution into a box integer plus a fractional part: x„ = Nn + 6n, 
where 

iVn+i = iV„ + a(^„) (15) 

On+l = f{On) (16) 

where a (6) = [f{6)] is the function returning the integer part of the argument. 

When considering the generating function (jHl) translation invariance of the system 
permits to restrict the average over the elementary cell. The transfer operator on the 
fundamental cell is written as: 

[C^pM = [ d^^eW»-^+'^W)5(/(^) _ 0)p(^^) (17) 



JT 

where all quantities refer to the torus map once we attach to every point in T the 
correspondent jumping number a{x) = [f{x)] jHllE]- 

The leading eigenvalue is the inverse of the smallest z solving the secular equation: 

det(l - zCp) = Fp{z) = (18) 

{Fj3{z) is usually called the spectral determinant). Using the formal relationship: 




det(l - zCfs) = exp Tr ln(l - zCp) = exp ( - — Tr Z:;^ ) (19) 
we are lead to computing the formal trace of C^: 

Tr£^ = Jdx 5{r{x) - a:)e(^"(")-"+""(")) = (20) 

„(/3o-„(x)) 

x=f"ix) 
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where 



n-l 



K{x)= n/'(/^(^)) (22) 



i=0 

n 

Y.^{rix)) (23) 

i=i 

and the sum runs over all periodic orbits of the system under study. 

We may now proceed to derive the standard expression for spectral determinants 
and zeta functions: the first step consists in using a geometric series expansion for the 
denominator in (PT|) : 

^'^'^r^ ^ |A„(x)|(A„(a;))^ ^^^^ 

fc=0 x=f"ix) 

We remark that the geometric series expansion makes sense only if we deal with a 
hyperbolic system, for which |A„(x)| > 1 for any point and iteration order (differently 
from what happens for intermittent systems). Next we perform a subtle step: the sum 
is over periodic points of any possible order n: so for instance fixed points appear at any 
order, and every p orbit reappears whenever p divides n. Moreover all points of a given 
periodic orbit contribute to the trace, but they all yield the same weight, as instabilities 
and jumping numbers cr„(x) are the same for every points of the same periodic orbit. 
So we may convert the sum over periodic points to a sum over orbits p, each taken once: 
we furher denote by Ap and dp respectively the instabilities and jumping numbers of the 
orbit p, each computed with respect to its prime period Up. In this way we get: 

{p} r=l k=0 ^ I^^pI J 

where {p} denotes the set of periodic orbits, each taken once, Up being the prime period 
of the orbit. The sum over r reproduces each repetition of a single orbit p in the original 
summation over periodic points: this sum may be performed exactly, leading to: 



det(l - zCf,) = exp V Vln 1 - — — = (26) 




fc=0 {p} ^ I PI p / 

Thus the spectral determinant may be expressed by an infinite product over the so called 
generalized zeta functions: 



{p} 



U(m (29) 



fe=0 
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and the leading zero can be calculated from the zero order zeta, the dynamical zeta 
function C(o)V 

We stress again that the quantities that enter the definition of zeta functions are 
the prime period Up of the orbit p, its instability Ap = YYiZo ^ f'iPi^p) ^^"^ integer 
factor (jp, that accounts for the orbit's behaviour once we unfold it on the real line. As 
a matter of fact, given any point Xp belonging to p {i.e. f"'^{xp) = Xp) we may either 
have that it is also a periodic point of the lift / {f^^{xp) = Xp), or it might be a running 
mode, i.e. f^^ixp) = Xp + ap, with ap G Z* (see (UlI)). 

It is well known that zeta functions calculation may be recast in terms of 
perturbative expansions (the order being determined by the longest prime period 
available), that yield exponentially good estimates In order to show that let's use a 
simplified example: suppose there is a complete grammar in the simple alphabet {0, 1}, 
calling for shortness tp = ^I'^'j^T we can expand the infinite product in series: 

{p} 

= (1 - to) ■ (1 - ti) ■ (1 - toi) ■ (1 - tool) ■ (1 - ton) ■ ■ ■ 
= 1 - to - ti + 

— (^01 — ^o^i) ~ (^001 ~ ^o^oi) ^ (^110 "~ ^i^oi) — ■ ■ ■ (31) 

We can notice that there are fundamental terms to and ti while others curvature 
term are made of a prime cycle shadowed by combinations of shorter cycles. 
Completely hyperbolic systems with a finite grammar yield curvature corrections that 
are exponentially small with the maximum period, so that ()31|) (power law expansion 
of the zeta function) is indeed a perturbative series. 

We now sketch the technique we will actually use to investigate arbitrary order 
transport moments: first we reorder the zeros Zi of the spectral determinants so that 
the dominant ones come first and write in general: 

F,{e-n = - ^^e-) (32) 



so that 



-^lnF^(e-^)= V-i^^^ (33) 



and if we now take the inverse Laplace transform, we get 

1 /-a+ioo 7 

''-"'^^'"[^'^('-'HE't (34) 

As the asymptotic behaviour is dominated by the leading eigenvalue we may extract 
it by using the dynamical zeta function instead of the spectral determinant, thus 
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Remember that by probability conservation we have -2o(0) = 1 and all transport 
properties are captured by small /3 expansion (fTT|) : 

gk 



(yk{n) = {{xn-xor)o = ^Gn{(3) 



i>a+ioo J 
J a-ioo ""^ '- 



13=0 

' ' (36) 



/3=0 



The evaluation of the integral on the right hand side of (j36|l requires dealing with high 
order derivatives of a composite function: this is generally accomplished by making use 
of Faa di Bruno formula: 

(7" ^ n ' d^ M 

-ff(LW)=i: j: ^^—(m).Bmt)) (37) 

B^[L{t))=[-—\ ■■■ 3-— (38) 



1! dt ) \n\ dt^ 
k = {ki, . . . kn} with ki = k, i ■ ki = n (39) 

Remember that the characterization of anomalous transport is captured by the 
function z/(g) (the fundamental quantity we want to evaluate) that expresses the 
asymptotic growth of moments of arbitrary order: 

{\xt~xon ~ r(^) (40) 

We already mentioned how ordinary diffusion yields h>{q) = q/2, while anomalous 
transport typically shows a non trivial behaviour, which cannot be encoded by a single 
exponent, but rather typically exhibits a phase transition [TH | I12 | ITT)]. In next sections we 
will apply the formalism to both normal and anomalous examples. These computations 
accomplish a twofold goal: besides checking the theory, they will also illustrate how 
subtle features of the underlying dynamical system are automatically included in the 
formalism. 



3. Normal transport 

We now show how the technique is able to reproduce the well known exponents q/2 for 
the transport moments of a completly chaotic dynamic. This oversimplified example is 
a warmup exercise before facing more involved (and interesting) cases. 
Consider the following map on the real line 

fix) = h'' ^ (41) 

^ \aix-l) + l a; G (1/2,1] 

For a generic value of the parameter a symbolic dynamics of pT|l is quite complicated 
(and this leads to interesting consequences, like fractal dependence of the diffusion 
constant on the slope a ^Zj): on the contrary by choosing an integer slope the torus 
map consists of complete branches: for example for a = 2m, m G N the corresponding 
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torus map consists of 2m complete branches and the grammar is thus unrestricted in a 
2m alphabet. Moreover the linearity of the map leads to complete cancellations of the 
curvature terms: whenever the weights associated to each cycle depend multiplicatively 
on the symbolic code entries, the contribution of a cycle and corresponding pseudocycles 
cancel out exactly and zeta functions pick up contribution only from fixed points. 

In this case there are 2m fixed points: the instabilities are trivially Aj = a for 
every fixed point while the jumping numbers associated to each branch are clearly: 
0,l,---,(a-l),-(a-l),---,-l,0. 

The zeta function for the torus map is then: 

QU^) = 1 - - ( 1 + E I (42) 

Long time behaviour of the q-th moment is captured by ()36|). Following jl] we 
concentrate first on (3 derivatives 



Faa di Bruno's expansion (jHTjl generates 



(43) 

/3=0 



k=i ki,-k„ ^om^^) 



(44) 

/3=0 

and a generic term contributing to the product in Bk{z) is now 

g^ji-W«('-)l.=o=|„»^" (45) 

Each term in ()44|1 has a different asymptotic behaviour and contribute differently to the 
transport moment once we perform the derivative with respect to s and take the inverse 
Laplace transform in ()36|) . The long time behaviour of the q-th transport moment is 
captured by choosing the leading singularity of the /3 derivatives near z = 1 (s = 0). 

The numerator in ()44|) doesn't diverge for z ^ 1 and the divergence comes only 
from the denominator: 

^ :i - z)-' (46) 



Co{/3}(^) 

The leading singularity comes from the term with the largest possible choice for k and 
hence all kj = 2 (remember that for kj odd the numerator vanish). This leads to the 
estimate 

gn 

and, performing the inverse Laplace transform, this yields z/(g) = q/2, typical of 
normally diffusing systems (where one is then interested in obtaining estimates of 
the prefactor, namely the diffusion constant). We turn now on more interesting case 
of intermittent diffusion where the presence of marginal points deeply modified the 
structure of transport exponents. 



lnCo-|mW|;3=o~(l-^)-"/' (47) 
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4. Intermittent systems 

In the last few years it lias been realized PHI ^1 1201 1^ 122] that weakly chaotic systems 
(in particular one dimensional intermittent maps, or infinite horizon Lorentz gas models 
lead to more complicated analytic structure of zeta functions, which typically exhibit 
branch points. In view of the inverse Laplace formula ()Hfi|l the modified analytic 
structure of the zeta function may induce anomalous behaviour (nonlinear diffusion 

[2311!). 

The qualitative difference between an hyperbolic system and an intermittent one 
can be well appreciated in fig. (0) where a long segment of an orbit for both cases is 
plotted against time. 




50 100 150 200 250 300 




50 100 150 200 250 300 



Figure 1. A segment of trajectory (300 iterations) witli the same initial condition 
is showed for two torus maps: a completely hyperbolic map (left) and one with 
intermittent fixed points (right). 

In the intermittent case chaotic wandering is interrupted by intervals of regular 
motion, generated by points that start their motion near the marginal point and slowly 
move away from it: we will see that there is a power low separation for nearby trajectories 
instead of the exponential growth we normally expect. In the next subsection we will use 
the thermodynamic formalism in order to detect the presence of power law behaviour 
in instabilities spectrum while in the following subsections we explain how it is possible 
to implement cycle expansion also in this anomalous case. 

4.1. Thermodynamic formalism 

In order to obtain information about the scaling behaviour of periodic orbits instabilities 
we turn our attention to a thermodynamic analysis of this set. Consider all periodic 
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points of n order i = 1 . . . N^. Consider now the partition of the unit interval with 
the intervals: 

Ii{n) = {Xi{n), Xi+l(n)) i = 1 . . . (48) 

whose length we denote by = |/i(n)|, and construct with them the following partition 
function: 

Zn{r) = J2 (kn)r^ (49) 

i=l 

where the parameter r is the formal analogue of an inverse temperature. By varying r 
one probes for all the relative different scaling of the intervals in the partition function 
sum For example for r — > +00 the partition function will be dominated by the 

thinnest intervals while for t —>■ ~oo the fattest ones prevail. 
Scaling properties may be analyzed by writing 

kn) = (iV„)-'^-(") (50) 

where we have introduced the local scaling indices it is important to notice that 

they will be finite if the interval shrink exponentially with n (as typically A^^ ~ a"') 
while will go to zero as (In n)/n for a power law behaviour. 
The (asymptotic) free energy is then defined as 

giT)= hm —L_\^Z^{t) (51) 

n^oo [n[JMn) 

ln(A^„) ^ 



n— »oo 



from which we can extract the average scaling index: 



Mr) = g\r) = hm ^^^^'X^ (53) 



Figure Q shows four approximated scaling index functions Unir) calculated for an 
intermittent map: this map has two complete branches and thus the number of periodic 
points of order n grows as 2" while the presence of marginal points generates a family of 
intervals that shrink polynomially with n (the exact form of this map will be presented 
as second example in following subsection (see ()7fj|l ). We can note that by increasing 
the hierarchical order the scaling index converge to finite values for r > while it 
slowly goes to zero for r < 0, as we expect due to intervals that shrink with power law 
behaviour. 

The thermodynamic analysis may be pushed further by rewriting the partition 
function as: 



Z„(r) = / d/iiV^^"(^)+'^") (54) 



mm 



where Nn"^^^^ is the density of intervals with scaling index ^. As the number of periodic 
orbits increases the integral can be approximated by the saddle point method: 

Nt'^' = C{n) \ - "(^)-'^) (55) 
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Figure 2. Approximated scaling index functions iJLn(j) versus r for the map in H76(l 
with 7 = 2.5 and with n=5,10,15,20. 



where is the solution of the extremum condition 
dSJr) 



d/i 



(56) 



and QniT) is the n-th order approximation to the free energy. 

The scahng spectrum s„(/i) is a highly irregular function of /i for finite n: in 
practise one evaluates Qnij) and replaces s(/i) = lim„_^oo ^^(/u) by its convex envelope 
^(/i), evaluated in the thermodynamic limit saddle point evaluation of (jSH). 

For the very same example used in calculating the scaling index of figure we now 
compute the convex envelope entropy S'„(/z). While for a complete hyperbolic dynamic 
we expect a convergence to a curve with strictly positive support, in our case marginality 
yields a tail in the small scaling region of S{^), which slowly advances towards zero as 
the hierarchical order increases, see figure Q: the right branch of entropy is stable 
and it is representative of exponential scalings while the leftmost branch slowly moves 
towards zero, being ruled by intervals shrinking with power laws. 



4-2. Transport 

We now apply the technique to two classes of one dimensional maps, where deviations 
from fully chaotic behaviour are provided by marginal fixed points of intermittent type 
[2n] • The prototype example of this class of systems is provided by the mapping on the 
unit interval 



where 7 > 1 is called the intermittency exponent. 



(57) 
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Figure 3. Approximated entropies 5'„(/i) versus /z for the map (|7?^ with 7 — 2.5 and 
with n=5, 10,15, 20. 

The presence of the marginal fixed point causes problems in using zeta function 
techniques in a straightforward manner: this fixed point has |A| = 1, which leads to 
divergences (see (12 The way out of this problem is to prune the fixed point away: 
from a symbolic dynamics point of view this leads to an infinite alphabet, since all 
0"1 fixed points lack curvature counterterms: the fundamental cycle part of the zeta 
function thus becomes (cfr. (jSH)) 

Cfund(^) = l-tl-toi- to^l to"! (58) 

SO even if curvature corrections are neglected we do not have a finite order polynomial: 
moreover the inclusion of curvature correction even in simple cases is by no means trivial 
[TH] . As we will see by considering specific examples, the analytic structure of will 
be strongly dependent on the intermittency exponent, through the behaviour of stability 
of cycles 0"1 coming closer and closer to the marginal fixed point. 

4-3. Intermittent diffusion I 

We can turn now on a specific example. The symmetry requirements are the same as 
the completely chaotic example (see relations (fT^ and (fT^ ). 

The map on the fundamental cell (see figure (j^) consists of two intermittent 
branches and one central hyperbolic region. If we assign jumping numbers ai = —1 to 
the left branch, crc = to the central one and = +1 to the right branch we get the 
lift map of figure (jlj). Whenever a particle remain trapped near a marginal point the 
corresponding unfolded trajectory on the real line consists of successive jumps on the 
neighbouring cell. 

The presence of complete branches leads to an unrestricted grammar in the symbolic 
code {L, C, R}, where the corresponding partition obviously consist of the inverse image 
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Figure 4. A torus map with intermittent fixed points, and the corresponding lift. 



of the unit interval under the corresponding three branches functions fi, fc^ fn- However 
we have to prune marginal fixed points L and R away, as explained before: the grammar 
is thus unrestricted on a countable new alphabet: 

{C, VC, UR, R^C, R^L} i, J, kJeN+ (59) 

The boundary points of such a partition can be constructed by sequences of iterates 
of the inverse map (L or R branches): /^^(l), /^^(l), . . . , /^^(O), /^^(O), . . .. The 
corresponding zeta function will thus consist of a (non polynomial) fundamental 
part ()58j] . together with curvature corrections: by a straightforward generalization of 
Gaspard-Wang piecewise linear approximation j2ZI of Pomeau-Manneville map, we can 
obtain an approximate model where all curvature terms vanish and the zeta function 
can be explicitly written without losing important dynamical features j2Hl- The power 
law scaling of the stability of cycles coming closer and closer to the marginal points can 
be extimated as 

A. ^ i-(»+^) (60) 

where i indicates the number of iteration the orbits needs to escape from the laminar 
phase, and 

« = 1/(7-1) (61) 

7 being the intermittency exponent of the L-R branches (see ()57|l ). 
The zeta function is: 

oo ^ 

C(o)y^) = l-az-bzY^ cosh(/?fc) (62) 

k=l 

where a and b are fixed by specifying the central region slope and the normalization 
condition C(o)o(-'-) ~ ^- 
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When equipped with information about the analytic structure of the zeta function, 
from (jH7|) we may single out the leading singularity in the logarithmic derivative of the 
zeta function, and then estimate the asymptotic behaviour of ()3fj|l . In our previous 
example the only divergence for 2; — >■ 1 in derivatives expansions ()37|) comes from the 
denominator while here the appearance of the Bose function (see (16 2j) ) 

9.i^) = E I (63) 
1=1 ' 

indeed may alter the analytic features. We recall how the behaviour as z 1~ depends 
upon fi (related to the intermittency exponent) 

(l-^)'^-^ fi < 1 

ln(l — z) fi = 1 

g,{z)^-{ C{f^)+C^{l-zr-' + D^{l-z) /xG(l,2) (64) 

C(2) + C2(l-2)ln(l-^) /i = 2 

afi) + C,{l-z) fi>2 

We stress again that this functions appear as a consequence of particular sequences 
of orbits coming closer and closer to the marginal fixed points whose stability increases 
only polynomially with the period [211111, a clear signature of local deviation from typical 
hyperbolic behaviour (which is ruled by exponential instability growth). 

Now look at a generic term in (ji^ and denote it by 2^fci...fc„- Taking into account 

that 



Q/qX(0)V^'' 



13=0 




' (65) 
i even 



we have that 



\^{IS)0\^)) j ^ki...k„ 

where the picks up the contributions from the product of Bose functions, and all j 
must be even, due to First we consider the case a e (0, 1), which corresponds to 

intermittency exponent 7 > 2: we have "D^^ ~ {1 — z)^°', that, together with 
implies that the dominant singularity is of the form 

~ (T^ (67) 
where p is determined by 

p = sup {ka + ^^(j — a)kj) = n (68) 

{ki...k„} j 

Once plugged into (jHUj) this leads to the estimate 

uiq) = q (69) 

which means that the whole set of transport moments is ruled by ballistic behaviour (at 
least for even exponents, where the method applies). We now turn to the more subtle 
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case a > 1: since the dynamical zeta function has a simple zero we get ~ (l — z)'', 

while the terms appearing in T>~^ modify the singular behaviour near z = 1 only for 
sufficiently high j 

If all {j} are less than a then the singularity is determined by V^: keeping in mind that 
the highest k value is achieved by choosing j = 2 and k2 = n/2, we get, by proceeding 
as before 

Kg) = 1 q<a (71) 

When q exceeds a we have to take into account possible additional singularities in 
and thus we get 

^" ~ why. (^2' 

where p is determined by 

/; ■ \, \ I "-72 n<2(a-l) 

p= sup A;+> {j-a)kj) = 'l ' o i 73 

{k,X} ^ ' \n+l-a n>2(a-l) 

which, once we take ()71|) into account, yields 

= I ' < t - '\ (74) 

[ g + 1 - a g > 2(a - 1) ^ ^ 

The set of exponents thus has a nontrivial structure, characterized by a sort of phase 
transition for q = 2{a — 1), a rather universal feature of many systems exhibiting 
anomalous transport jl2^ We notice that the parameter ruling the presence of a 
phase transition (and the explicit form of the spectrum) is a, that is the exponent 
describing the polynomial instability growth of the family of periodic orbits coming 
closer and closer to the marginal fixed point, and thus describing the sticking to the 
regular part of the phase space: in the present example a = (7 — 1)"^, and thus 
the sticking exponent is easily connected to the intermittency index. We remark that 
our reference to "phase transitions" is somehow arbitrary, since in effect our technique 
applies in the present form only to even moments: however results as ()74|1 are consistent 
with all our numerical checks for the whole real q axis. 



4.4- Intermittent diffusion II 

We now discuss a further example, which will provide additional evidence of the virtues 
of our technique. While the setting is similar to the former case we slightly modify it to 
conform to original references. The map we consider belongs to a class of 2-dimensional 
systems were the diffusive process takes place in a direction while all the dynamics is 
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limited to the torus on the other dimension: 

Xt+l = f{xt) 

Vt+i = yt + xt 



(75) 



where f{x) is a torus map, again dependent on an intermittency parameter 7, imphcitly 
defined on T = [—1, 1) in the following way [TH] : 



X 



2^ 



fix) + ^ ( 1 - fix) 



0<x < 1/(27) 
1/(27) < a; < 1 



(76) 



for negative values of x the map is defined as /(— x) = —fix) (see (O): the map 
is plotted in fig. Diffusion can be treated by the same techniques of the former 

subsection if we are able to control the values the function (equivalent to the jumping 
numbers of previous example) 



E 

t=i 



Xf 



(77) 



takes on the periodic orbits. 

The most striking difference with respect to the former example is that, for any 
values of J, the map (ffF)|) has a constant invariant measure [TK] . 



f(x) 




Figure 5. The intermittent map witii constant invariant measure. 

This comes from the summation property 
1 



y=f-Hx) 



f'iy) 



1 



in 
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that easily follows once we observe that each point has exactly two counter images 
yi{x) and y2{x): they come from two different pieces in ()76|) and they have opposite 
sign: by deriving the two pieces of ()7fi|l then we get the summation rule. This is quite 
different from the former case where the torus map has non trivial ergodic properties: an 
absolutely continuous invariant measure only exists for a > 1 (see for instance jSHI and 
references therein) while the ergodic behaviour is much more complex when a < 1. In 
any case sticking induces peaking of the measure around marginal fixed points (see \27\). 
In this sense the behaviour of ()7fi|l is similar to what happens in an area preserving map 
example [3U], where a parabolic fixed point coexist with a Lebesgue invariant measure. 

The difference with the previous example can be euristically understood if we note 
that now the slope of the map in the chaotic region is not bounded from above: the 
closer we need to be injected near the marginal points the lesser the probability that this 
happens. This feature deeply modifies the instabilities of periodic orbits coming closer 
and closer to the marginal fixed points and we expect that a linearization as before will 
not be able to reproduce them. 

A linearization is still possible but we must put attention on matching the condition 
()78|) . To do that let's start by identifying for each branch a laminar and an injection 
regions (for example we call B = [0, 1/ (27)] and B = [1/ (27), 1] respectively the injection 
and laminar regions of the right branch (see figure ©)• 

A, A^.. B, B_,Bj 



























1 

1 








1 






i 


1 



A3 A, A, °..B,B, 



Figure 6. Linearization of the map (|76|l . Near marginal points, laminar regions, a 
Gaspard-Wang type partition is used while a finer one is used in the injection zone . 

The Gaspard-Wang linearization has to be modified in the injection zones A and 
B: a linear approximation for all such intervals will give an uniform probability of 
being reinjected in the laminar zone and this does not respect the peculiar ergodic 
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property of this map. The way out to this problem is to perform a finer partition 
of that regions taking the inverse image (inside the injection zones) of each interval 
partitioning the corresponding laminar zone. For example we will call Bi the injection 
interval corresponding to inverse image of the laminar one Ai, B2 corresponding to A2 
and so on (see figure ©)• The linearization is now straightforward for each of this 
intervals. Moreover the instabilities of all the periodic orbit are well reproduced once 
we note that the slope of the linearized map restricted to such injection intervals can be 
safely approximated by (see summation property (|7Hjl ): 

= 1 - Aa! (79) 

We will now show how such a linearization is able to reproduce the behaviour of 
periodic orbit instabilities. We have seen that the periodic orbits that are relevants for 
the construction of the zeta function are the ones that accumulate near the marginal 
fixed point and for which the instability has a power law behaviour. In this new map 
the instability of such orbits take a linear correction factor (with respect to the length) 
from the injection zone: the intervals in the laminar regions, constructed in the same 
spirit of Caspar d- Wang linearization [2Zj, shrink polinomially with the index k: 

I'^-j^n (80) 
from wich is simple to derive the periodic orbit instabilities. We rewrite the instabilities 
of the relevant orbits for the previous example (I) in comparison with the ones of this 
map (II): 

Al.„ = lj^^ji...^l= (81) 

A;"+^ (82) 



A^^.. = [|^|i...^][l-y^r[l-|ir (83) 
-A[.^-^^~A:Ai.^ (84) 

where we note that for this second example the modified linearization produce two new 
factors (see (jEHI)) that correspond to the slope of the two injection intervals visited by 
the orbit; one of them depends explicitly on the lenght of the orbit and gives a relevant 
correction to the instability of the cycle (j84p . 

This trend is what we expect from numerical simulations (see figure (|7j)). Note that 
this new linearization is able to capture the behaviour of the instabilities of all orbits, 
not only the one accumulating to the marginal point. As a further example we may 
consider periodic orbits of the type A^B^ were the power law exponent is twice the one 
before, in agree with performed numerical simulations. Thus fl78|l deeply modifies the 
relationship between the intermittency exponent and the instabilities of periodic orbits 
shadowing the marginal fixed points. 
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3 - 




4 



a 



Figure 7. The exponent x(q^) of power law decay of instabilities ^^Jkg ^ fc-^'"-* as a 
function of the parameter a (open circles) and the best fit xi'^) = (0.88a + 1.9) (full 
line). 



Once we notice that, for the cycles dominating the zeta function expression, the 
jumping function ()77j) can be approximate by the lenght of the periodic orbit (with the 
appropriate sign), the piecewise hnear approximation yields the dynamical zeta function 



where again a = 1/(7 — 1). 

By repeating the steps of the calculation performed in the former case we get that 
we always have a phase transition with 



which may also be checked numerically (see figure (|H))). 
4-5. The Lorentz gas with infinite horizon 

The technique does not only apply to one dimensional systems, but is capable of dealing 
with more complicated systems, if we are able to characterize the sequence of orbits 
that probe more and more accurately the local behaviour near the ordered region. A 
remarkable example in this respect is the Lorentz gas with infinite horizon (see fig. (jH))), 
where a particle collides elastically with a regular array of fixed circular scatterers [3], 
the corresponding torus map being a square billiard table with a circular scatterer in the 
center (the so called Sinai billiard In this example the dynamical feature that 

lowers the chaoticity of the system is provided by the possibility of infinite free flights 
along corridors (for instance horizontal and vertical direction in fig. (jHl)). Periodic 
orbits for the torus maps (running modes for the full system) that mimic the free flight 




(86) 




q <2a 
q> 2a 



(87) 
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Figure 8. Spectrum of the transport moments for the map l|76() with 7 = 1.5: the best 
fit on numerical data is ^{q) = 0.50q + 0.04 for the dotted line, and i'{q) = 0.98(7 — 1-82 
for the full line. 



are for instance orbits that bounce from the (m, 0) row to the {n, 1) row travelhng a 
distance of p lattice spacings between each two successive colhsions. In view of our 
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Figure 9. Unbounded horizon Lorentz gas, with discs labelling. 

former considerations the estimate we need is how the stabihty of such orbits scales in 
the big p limit. Such a calculation is equivalent to computing the fraction of orbits that 
starting from the (0, 0) disc suffer their first collision exactly in the {p, 1) disk: such a 
probability decays like 1/p^ (see [20]), and in this sense the Lorentz gas corresponds to 
the case a = 2, in This leads to a direct estimate of the spectrum of exponents, 

according to ((ZH): 

It is interesting to show how the phase transition sits exactly at g = 2, that is the order 
corresponding to diffusion. The behaviour described in ()88|1 . was obtained recently in 
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[TB] (actually for a wider class of extended billiards), by different methods. We also 
remark that our method is also capable of taking track of logarithmic corrections to 
the asymptotic estimates, and in the present example we get that the variance (second 
moment) asymptotically grows as t ■ Int (see [TIT]). The appearance of a logarithmic 
correction may be understood if we notice that, for a = 2, we have that 

5^W^ )~ (l-e-) 

(cfr. ()(i2j) ): then the leading singularity, dominating the integral is of the form 
ln(l — e~'*)/(l — e~*)^, which, by using for instance Tauberian theorems for Laplace 
transforms, leads to the logarithmic correction in time. Notice that our results are 
asymptotic and they do not depend on the geometric features of the billiard (radius to 
cell size ratio), which would of course modify transient behaviour, as well as prefactors 
(with dramatic effects once this ratio tends to or 1). 

5. Conclusions 

In this paper we have proposed a zeta function technique to compute the asymptotic 
behaviour of different transport moments: the essential ingredient in the analysis of 
weakly chaotic examples, that show deviations from the normal diffusing case, has been 
shown to be a proper characterization of the sequence of periodic orbits probing closer 
and closer dynamical features of the marginal structures. In particular the behaviour of 
the moments' spectrum z/(g) is governed by the power law relating such orbits' instability 
to their period. In many cases i/(g) is characterized by a phase transition between a 
normal regime and a ballistic one. 
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